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Abstract 

The computation of the correction Ar in the W-Z mass correlation, derived from muon 
decay, is described at the two-loop level in the Standard Model. Technical aspects which 
become relevant at this level are studied, e.g. gauge-parameter independent mass renormal- 
ization, ghost-sector renormalization and the treatment of 75. Exact results for Ar and the W 
mass prediction including 0{a^) corrections with fermion loops are presented and compared 
with previous results of a next-to-leading order expansion in the top-quark mass. 



*Talk given by A. F. at the 5th Zeuthen Workshop on Elementary Particle Theory "Loops and Legs in Quantum 
Field Theory", Bastei/Konigstein, Germany, April 9-14, 2000 



1 Introduction 



The electroweak Standard Model (SM), together with the theory of strong interactions (Quantum 
Chromodynamics, QCD), provides a comprehensive description of experimental data with remark- 
able consistency. In order to further test the validity of the SM predictions and restrict its only 
missing parameter, the Higgs boson mass, precision measurements play a key role. The analysis 
of precision observables is sensitive to quantum corrections in the theoretical predictions, which 
depend on all the parameters of the model. In this way, the top-quark mass, nit, had been predicted 
in the region where it was experimentally found. 

The constraints on the Higgs mass, Mh, are still rather weak since Mh appears only loga- 
rithmically in the leading order SM predictions. Therefore it is of high interest to further reduce 
experimental and theoretical uncertainties. An important quantity in this context is the quantum 
correction Ar in the relation of the gauge boson masses M^, Mz with the Fermi constant Gp and 
the fine structure constant a [I]. 

This relation is established in terms of the muon decay width, which was first described in the 
Fermi Model as a four-fermion interaction — > e" i/^ with coupling Gp. This yields for the 
muon decay width 

where F{ml/mj^) subsumes effects of the electron mass on the final-state phase space and Aq de- 
notes the QED corrections in the Fermi Model. Including two-loop QED corrections [2, 3] one ob- 
tains from the measurement of the muon decay width [4] Gp = (1.16637 ± 0.00001) 10~^ GeV~^. 

The decay process in the SM, on the other side, involves the exchange of a W boson and 
additional electroweak higher order corrections. Although not part of the Fermi Model, tree-level 
W propagator effects are conventionally included in eq. (1) by means of a factor (1 + |m^/M^). 
Yet, this is of no numerical significance. The relation between the Fermi constant Gp and the SM 
parameters is expressed as 

Here the SM radiative corrections are included in Ar (with = 1 — M^/M^). 

Since Gp is known with high accuracy it can be taken as an input parameter in equation (2) in 
order to obtain a prediction for the W mass (M^ is still afflicted with a considerable experimental 
error, Mw = 80.419 ± 0.038 GeV [5]): 




1 + Ar 



V2GpMl 



(3) 



Since Ar itself depends on Mw, eq. (3) is to be understood as an implicit expression. 

One expects substantial improvement on the W mass determination from future colliders, 
which will allow for increasingly stringent constraints on the SM from the comparison of the pre- 
diction with the experimental value for Mw. This requires an accurate theoretical determination of 
the quantity Ar. 



1 



The one-loop result [1,6] can be split into the following contributions: 

Ar(") = Aa - + ^r,,^{M^). (4) 

Dominant corrections arise from the shift in the fine-structure constant, Aa, due to large logarithms 
of light-fermion masses i;^ 6%), and from the leading contribution oc to the p parameter result- 
ing from the top/bottom doublet, which enters through Ap (a; 3.3%). The full Mn-dependence is 
contained in the remainder Arrem (~ !%)• 

Furthermore, QCD corrections of 0{aa^) [7] and 0{aa1) [8] are known. Leading electroweak 
fermionic 0{a'^) contributions were first taken into account by means of resummation relations [9]. 
Different approaches for the calculation of electroweak two-loop corrections have been pursued, 
involving expansions for large values of Mh [10] and mj [11, 12, 13]. 

It turned out that both the leading 0{a'^mf/M^) [12] and the next-to-leading 0{a'^m'^/M^) 
[13] coefficients in the mt expansion yield important corrections of comparable size. Therefore a 
complete two-loop calculation of fermionic contributions would be desirable in order to further re- 
duce the theoretical uncertainty, in particular if one considers that these contributions are dominant 
already in the one-loop result. 

A first step into this direction was the determination of the exact Higgs mass dependence of the 
fermionic 0{a'^) corrections to Ar [15]. Recently the full calculation of these contributions [16] 
has been accomplished. The results include all electroweak two-loop diagrams with one or two 
fermion loops without any expansion in mt or Mh. This talk gives an overview on the techniques 
and the results of this calculation. 



2 Calculational methods 

Since in the present calculation all possibly infrared divergent photonic corrections are already 
contained in the definition (1) of the Fermi constant Gp and mass singularities are absorbed in the 
running of the electromagnetic coupling, Mw represents the scale for the electroweak corrections 
in Ar. Therefore it is possible to neglect all fermion masses except the top quark mass and the 
momenta of the external leptons so that the muon decay diagrams reduce to vacuum diagrams. 

All QED contributions to the Fermi Model have to be excluded in the computation of Ar as 
they are separated off in the definition of Gf, see eq. (1). When writing the factor (l+Aq) formally 
as (1 -I- Ao;)^ the decay width up to 0{a'^) can be decomposed as 



= r^,tree(l + Aujf (1 + A? 



^2 



F 

/i,tree 



1 + 2(Acj(«) + Ar(")) + (Acu(") + Ar^"))^ + 2Aa;(")Ar(") (5) 
+2Ac^("') + 2Ar("') + 0{a^ 



Apart from the one-loop contributions this includes two-loop QED corrections Acu^"'^^ and mixed 
contributions of QED and weak corrections Acu*^"^ Ar*^"^ which both thus have to be excluded 
in Ar'^" \ For fermionic two-loop diagrams it is possible to find a one-to-one correspondence 
between QED graphs in Fermi Model and SM contributions. 
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The renormalization is performed in the on- shell scheme. In this context the mass renormal- 
ization constants require the computation of two-loop self-energy diagrams with non-vanishing 
momentum. 

All decay amplitudes and counterterm contributions have been generated with the program 
FeynArts 2.2 [17]. The amplitudes are algebraically reduced by means of a general tensor inte- 
gral decomposition for two-loop two-point functions with the program TwoCalc [18], leading to a 
fixed set of standard scalar integrals. Analytical formulae are known for the scalar one-loop [19] 
and two-loop [20] vacuum integrals whereas the two-loop self-energy diagrams can be evaluated 
numerically by means of one-dimensional integral representations [21]. 

In order to apply an additional check the calculations were performed within a covariant 
gauge, which introduces one gauge parameter i = 7,Z,W, for each gauge boson. It has been 
explicitly checked at the algebraic level that the gauge parameter dependence of the final result 
drops out. 



3 On-shell renormalization 



For the determination of the one-loop counterterms (CTs) and renormalization constants the con- 
ventions of ref. [22] are adopted. Two-loop renormalization constants enter via the counterterms 
for the transverse W propagator and the charged current vertex: 



w w 

.AAAA((WV^^AAA.<WW^ 



(6) 




V2s 



sz, 



Ss. 



e(2) 



W(2) 



+ 



{6Zl 



eh 
(2) 



(7) 



-I- (1-loop renormalization constants) 



Here, 6Z'^ denotes the field-renormalization constant of the field 0, SM^ the corresponding mass 
CT, and dZe the charge-renormalization constant. The numbers in parentheses indicate the loop 
order. The mixing angle CT, (5sw(2), is expressible through the gauge boson mass CTs. Throughout 
this paper, the two-loop contributions always include the subloop renormalization. 

The on-shell masses are defined as the position of the propagator poles. Starting at the two- 
loop level, it has to be taken into account that there is a difference between the definition of the 
mass as the real pole of the propagator D, 



and the real part M of the complex pole. 



0, 



= M -iMT. 



(8) 



(9) 



The imaginary part of the complex pole is associated with the width T. The definition (9) yields 
for the W mass CT 



(10) 
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whereas for the real pole definition the last term of eq. (10) is missing. E^, denote the 
transverse W self-energy and its momentum derivative. Similar expressions hold for the Z boson. 

The W and Z mass CTs determine the two-loop mixing angle CT, 5sw(2), which has to be 
gauge invariant since Sw is an observable quantity. With the use of a general gauge it could be 
explicitly checked that 5sw(2) is gauge-parameter independent for the complex-pole mass definition, 
whereas the real-pole definition leads to a gauge dependent 6s^^2)- This is in accordance with the 
expectation from S-matrix theory [23], where the complex pole represents a gauge-invariant mass 
definition. 

It should be noted that the mass definition via the complex pole corresponds to a Breit-Wigner 

parameterization of the resonance shape with a constant width. For the experimental determination 
of the gauge boson masses, however, a Breit-Wigner ansatz with a running width is used. This has 
to be accounted for by a shift of the values for the complex pole masses [24], 

_ 

M = M , (11) 

2M' ^ ^ 

which yields the relations 

Mz = Mz-34.1MeV, 

Mw = Mw - 27.4 (27.0) MeV for Mw = 80.4 (80.2) GeV. ^ ^ 

For Mz and Fz the experimental numbers are taken. The W mass is a calculated quantity, and 
therefore also a theoretical value for the W boson width should be applied here.' The results above 
are obtained from the approximate, but sufficiently accurate expression for the W width. 



\ = 1 + . (13) 



At the subloop level, also the Faddeev-Popov ghost sector has to be renormalized. The gauge- 
fixing sector for the gauge fields A^, Z^, W^^ (x, 0=*= denote the unphysical Higgs scalars) 



= - 1 ((F^)2 + (F^)2 + F+F- + F-F+) , 



F^^iay'-d^A^ + ^d^Z^, (14) 



= (ef ) d,Z^ + ^d,A>^ - (ef ) ^ Mz X, 

does not need renormalization. Accordingly, one can either introduce the gauge-fixing term af- 
ter renormalization or renormalize the gauge parameters in such a way that they compensate the 
renormalization of the fields and masses. Both methods ensure that no counterterms arise from the 

^In the version of this paper appearing in the proceedings for simplicity a fixed shift of leading order has been used. 
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gauge-fixing sector but they differ in the treatment of the ghost Lagrangian, which is given by the 
variation of the functionals F" under infinitesimal gauge transformations 56'^, 

S (15) 

o,6=7,Z,± 

In the latter case, which was applied in this work, additional counterterm contributions for the ghost 

sector arise from the gauge parameter renormalization. The parameters in (14) are renormalized 
such that their CTs 5^° exactly cancel the contributions from the renormalization of the fields and 
masses and that the renormalized gauge parameters comply with the gauge. 

4 Treatment of the 75-problem 

In four dimensions the algebra of the 75-matrix is defined by the two relations 

{75,7a} = for a = l,...,4 (16) 
Tr{757'^7"7''7''} = Ue'"''"'. (17) 

It is impossible to translate both relations simultaneously into D ^ A dimensions without encoun- 
tering inconsistencies [25]. 

A certain treatment of 75 might break symmetries, i.e. violate Slavnov-Taylor (ST) identities 
which would have to be restored with extra counterterms. Even after this procedure a residual 
scheme dependence can persist which is associated with e-tensor expressions originating from 
the treatment of (17). Such expressions cannot be canceled by counterterms. If they broke ST 
identities this would give rise to anomalies. 

't Hooft and Veltman [25] suggested a consistent scheme which was formulated by Breit- 
enlohner and Maison [26] as a separation of the first four and the remaining dimensions of the 
7-Matrices (HVBM- scheme). It has been shown [27] that the SM with HVBM regularization is 
anomaly-free and renormalizable. This shows that e-tensor terms do not get merged with diver- 
gences. 

The naively anti-commuting scheme, which is widely used for one-loop calculations, extends 
the rule (16) to D dimensions but abandons (17), 

{75,7a} = for a^l,...,D (18) 

Tr{757^7'^'"7"} = 0. (19) 

This scheme is unambiguous but does not reproduce the four-dimensional case. 

In the SM particularly triangle diagrams (Fig. 1) containing chiral couplings are sensitive to 
the 75-problem. For the presented work one-loop triangle diagrams have been explicitly calcu- 
lated in both schemes. While the naive scheme immediately respects all ST identities the HVBM 
scheme requires the introduction of additional finite counterterms. Even after this procedure fi- 
nite differences remain between the results of the two schemes, showing that the naive scheme is 
inapplicable in this case. 
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Figure 1: Generic Standard Model triangle diagrams subject to the 75-problem. 




Figure 2: CC vertex diagrams with triangle subgraphs. 



In the calculation of Ar triangle diagrams appear as subloops of two-loop charged current 
(CC) vertex diagrams (Fig. 2). One finds that for the difference terms between both schemes this 
loop can be evaluated in four dimensions without further difficulties. This can be explained by the 
fact that renormalizability forbids divergent contributions to e-tensor terms from higher loops in 
the HVBM scheme. The e-tensor contributions from the triangle subgraph in the HVBM scheme 
meet a second e-tensor term from the outer fermion lines in Fig. 2, thereby resulting in a non-zero 
contribution to Ar. 

Computations in the HVBM scheme can get very tedious because of the necessity of additional 
counterterms. Therefore another method shall be examined. One can consider a "mixed" scheme 
that uses both relations (16) and (17) in dimensions despite their mathematical inconsistency. 
This scheme is plagued by ambiguities of 0(D — 4). When applied to the calculation of one- 
loop triangle diagrams the results immediately respect all ST identities and differ from the HVBM 
results only by terms of 0{D — A), 

r^ir = rs'fi) + op-4). (20) 

Since for the difference terms the second loop can be evaluated in four dimensions, this also holds 
for the two-loop CC diagrams, 

r?c(T = r-g2) + (^(^-4). (21) 

Thus the mixed scheme, despite being mathematically inconsistent, can serve as a technically 
easy prescription for the correct calculation of the CC two-loop contributions. 

5 Results 

In the previous sections the characteristics of the calculation of electroweak two-loop contributions 
to Ar have been pointed out. Combining the fermionic 0{a'^) contributions with the one-loop and 
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Figure 3: Contribution of one-loop and higher order corrections to Ar. 



the QCD corrections yields the total result 

Ar = Ar^") + Ar^""^) + Ar^"""') + Ar^^^^'^ + Ar^^/""'^. (22) 

Here A^j, Nj symbolize one and two fermionic loops respectively. Fig. 3 shows that both the QCD 
and electro weak two-loop corrections give sizeable contributions of 10-15% with respect to the 
one-loop result. 

In Fig. 4 the prediction for derived from the result (22) and the relation (3) is compared 
with the experimental value for Mw. Dotted lines indicate one standard deviation bounds. The 
main uncertainties of the prediction originate from the experimental errors of rrit — (174.3 ± 5.1) 
GeV [4] and Aa = 0.05954 ± 0.00065 [28]. It can be noted that light Higgs masses are favored 
by this analysis. Further implications of the precision calculation of A/w are discussed in [30]. 

These results can be compared with the results obtained by expansion of the two-loop contri- 
butions up to next-to-leading order in mt [13, 14]. The predicted values for for several values 
of Mh are given in Tab. 1. Agreement is found between the results with maximal deviations of 
about 4 MeV in Mw. 

The theoretical uncertainty due to missing higher order contributions can be estimated as fol- 
lows. The missing 0{a'^) purely bosonic corrections can be judged by means of resummation 
relations to be very small (< 1 MeV effect on M^ for a light Higgs). An estimate of the 0{a^) 
terms can be obtained from the renormalization scheme dependence of the two-loop result yield- 
ing about 2-3 MeV, and the missing higher order QCD corrections were estimated to be about 4-5 
MeV [31, 32]. Adding this up linearly, one arrives at a total uncertainty of about 7 MeV for the 
Mw-prediction at low Higgs masses (Mh< 150 GeV). 
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Figure 4: Prediction for Mw as function of Mr compared with the experimental W mass. 

6 Conclusion 



In this talk the realization of an exact two-loop calculation of fermionic contributions in the full 
electroweak SM and its application to the precise computation of Ar were presented. Numerical 
results and an estimate of the remaining uncertainties were given and might serve as ingredient for 
future SM fits. 



The speaker would like to thank S. Bauberger and D. Stockinger for valuable discussions. 



Table 1: 


Comparison between Afw- 


-predictions from a NLO expansion 


in nil {M^^'^) and the full 


calculation (M^"). dM^ denotes the difference. Experimental input values are taken from [14]. 


Mh 


M^^P^ 


M^" 


5Mw 


[GeV] 


[GeV] 


[GeV] 


[MeV] 


65 


80.4039 


80.3997 


4.2 


100 


80.3805 


80.3771 


3.4 


300 


80.3061 


80.3051 


1.0 


600 


80.2521 


80.2521 


0.0 


1000 


80.2129 


80.2134 


-0.5 
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